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ABSTRACT 

Context. The chromosphere of the Sun is a temporally and spatially very varying medium for which the assumption of ionisation 
equilibrium is questionable. 

Aims. Our aim is to determine the dominant processes and timescales for the ionisation equilibrium of calcium under solar chromo- 
spheric conditions. 

Methods. The study is based on numerical simulations with the RADYN code, which combines hydrodynamics with a detailed so- 
lution of the radiative transfer equation. The calculations include a detailed non-equilibrium treatment of hydrogen, calcium, and 
helium. Next to an hour long simulation sequence, additional simulations are produced, for which the stratification is slightly per- 
turbed so that a ionisation relaxation timescale can be determined. The simulations are characterised by upwards propagating shock 
waves, which cause strong temperature fluctuations and variations of the (non-equilibrium) ionisation degree of calcium. 
Results. The passage of a hot shock front leads to a strong net ionisation of Ca 11, rapidly followed by net recombination. The relax- 
ation timescale of the calcium ionisation state is found to be of the order of a few seconds at the top of the photosphere and 10 to 30 s 
in the upper chromosphere. At heights around 1 Mm, we find typical values around 60 s and in extreme cases up to ~ 150 s. Generally, 
the timescales are significantly reduced in the wakes of ubiquitous hot shock fronts. The timescales can be reliably determined from 
a simple analysis of the eigenvalues of the transition rate matrix. The timescales are dominated by the radiative recombination from 
Calll into the metastable Call energy levels of the 4d^D term. These transitions depend strongly on the density of free electrons and 
therefore on the (non-equilibrium) ionisation degree of hydrogen, which is the main electron donor. 

Conclusions. The ionisation/recombination timescales derived here are too long for the assumption of an instantaneous ionisation 
equilibrium to be valid and, on the other hand, are not long enough to warrant an assumption of a constant ionisation fraction. 
Fortunately, the ionisation degree of Call remains small in the height range, where the cores of the H, K, and the infrared triplet lines 
are formed. We conclude that the difference due to a detailed treatment of Ca ionisation has only negligible impact on the modelling 
of spectral lines of Ca II and the plasma properties under the conditions in the quiet solar chromosphere. 

Key words, hydrodynamics; shock waves; Sun: chromosphere; waves; Radiative transfer 
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1. Introduction rate equations instead of the much simpler statistical equilib- 

rium equations. In paper I, it was shown that the resulting 

^ Modern observations unambiguously show that the chromo- ^^^^^^^^ hydrogen ionisation degree is significantly re- 

^. sphere of the Sun is very dynamic and inhomogeneous (see ^^^^^ ^j^^ recombination rates. The ionisation degree 

the reviews by, e.g., |S .c hriiver 20^ Uudge 2006; Rutten 2007; ^j^^^^ ^^^^ j^^^j temperature, density, 

I Wedemeyer-Bohmetal. , ,2(^09.) The interpretation of chromo- ^j^^ ^.^^i^^i^j^ ^^j^ ^j^^ 

on the history of these proper- 

spheric observations which often involves the construction of ^.^^ ^j^.^^^ Rammacher & Cuntz ( 1991) considered the non- 
numerical models, thus must take into account spatial and tem- equilibrium ionisation of magnesium (Mg) and calcium (Ca) in 
poral variations. Unfortunately many simplifying equilibrium ^^^-^ ^^j^^ chromosphere models. Their numerical im- 
assumptions that can be made for the low photosphere, are not igj^entation included a non-local thermodynamic equilibrium 
applicable for the layers above. A realistic chromosphere model ^^lTE) approximation with very simplified two-level atoms for 
^ must therefore account for deviations from the equilibrium state jj jjj jj jjj 

- - in the context of an atmosphere that changes on short timescales t . • • ■ , ^ tt ttt • • ■ , , 

■ 11 1 1 A f 1 1 ■ »u 1 In this paper, we investigate the Ca II-III ionisation balance 

and small spatial scales. A successful example is the explana- . , ., rr,, . • , , .t^-i^xt • , • • 

c c 1 • • 5> u IT^i — i D ■ I in detail. The study IS based on RADYN simulations similar to 

tion or the formation of calcium grains by ICarlsson & Steinl , , , . t rr,, • , . , . 

jinmK f^^u 1 f • 1 ,1; J J , the hydrogen case in Paper I. The numerical simulations are 

(II997I) . Other examples of simulations with a time-dependent , , . „ Fn r,,, , r , ■ 

■ i-u ■ t * t • 1 J »u • • e u J described in Sect. Kl The results of this study are presented in 

non-equmbrium treatment include the ionisation of hydro- _ m , , .• , , • • o m 

,17!; — p o> ■ ' 'tTatio u f J . n T Sect. Qi followed by a discussion and conclusions in Sect. MJ 

gen (Carlsson & Stein 2002, hereafter referred to as Paper I; ' -' ' 

iLeenaarts & Wedemeyer-Bohm 2006; Leen aarts et al. 20 07) and 

the concentration of carbon monox ide dAsen sio Ram os et al.l 

l2003t IWedemever-Bohm et al.ll2005h . When an individual pro- 2. Numerical simulations 

cess, be it a chemical reaction or atomic transition, changes 

on finite timescales longer than the dynamic timescales of the We use one-dimensional radiation hydrodynamic simulations 

atmosphere, then it becomes necessary to solve a system of calculated with RADYN, which ar e in most aspect s very similai - 

to the earlier simulation runs by Carlsson & Steinl d 19921 11994 

Send ojfprint requests to: sven.wedemeyer-bohm@astro.uio.no I1995L paper I). Here, we summarise only the most important 
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Fig. 1. Calcium model atom with the five lowest energy levels 
of Ca II, a continuum level (Ca III), and all included b-b and b-f 
transitions. Energies are not to scale. 



properties. More details can be found in paper I and references 
therein. 

The simulation code RADYN solves the equations of mass, 
momentum, energy, and charge conservation together with the 
non-LTE radiative transfer and population rate equations, on 
an adaptive mesh in one spatial dimension. It takes into ac- 
count non-equilibrium ionisation, excitation, and radiative en- 
ergy exchange from the atomic species H, He, and Ca with 
back-coupling on the hydrodynamics. Also, the effect of mo- 
tion on the emitted radiation from these species is considered. 
Hydrogen and singly ionised calcium are modelled with six- 
level atoms and helium with a nine-level atom. Furthermore, 
doubly ionised helium is included. The transitions between all 
considered atomic levels are treated in detail. Each line is de- 
scribed with 31-101 frequency points, whereas 4-23 frequency 
points are used for each continuum. Other elements than H, He, 
and Ca are taken into account in the form of background con- 
tinua in LTE, which are de rived with the Uppsala atmospheres 
program ( lGustafssonlll973h . 

The boundary conditions are the same as for the simulations 
in paper I. The lower and upper boundaries are both transmit- 
ting. The lower boundary is located at a fixed geometrical depth, 
which corresponds to 480 km below T500 = 1 in the initial at- 
mosphere. We again use the same piston as in paper I to excite 
waves at the bottom, which then propagate through the atmo- 
sphere. The piston velocity is based on Doppler-shift observa- 
tions in a Fe I line at A = 396.68 nm in the wing of the Ca H line 
dLites et al.lll993h . The upper boundary condition is located at a 
height of 10 Mm. It represents a corona on top of the simulated 
layers with a te mperature set to 10'' K and corresponding inci- 
dent radiation (lTobiskalll99ll) . The dynamic simulation has an 
overall duration of 3600 s. 

The numerical simulation is characterised by shock waves, 
which are excited as acoustic disturbances of small amplitude 
through the piston at the lower boundary. They propagate up- 
wards and steepen into shocks with a saw-tooth profile above the 
photosphere. The simulated chromosphere is therefore subject to 
strong fluctuations in all hydrodynamic quantities. It takes a few 
100 s until the first shock waves transform the initially hydro- 
static atmosphere into the characteristic dynamic atmosphere. 
The first 600 s of the simulation are therefore excluded from the 
analysis. 
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Fig. 2. Ionisation fractions of calcium (xcn, thick solid, see 
Eq. ([T]i) and of hydrogen (dashed) as a function of column mass 
in the initial radiative equilibrium atmosphere. The height scale 
Z() of the initial atmosphere is given as reference. The Call-III 
ionisation fraction rises from a minimum of 2 lO""* at Ig - 0.0 
(z ~ 180 km) in the low photosphere to ~ 5% at the base of 
the transition region. The grey lines represent ;t'ca in simulation 
snapshots at an interval of 10 s during the whole simulation se- 
quence of 3600 s. The thin short-dashed lines are the 5 % and 
95 % percentiles. 



A large number of additional short simulation runs are car- 
ried out, which are used for a timescale analysis (see Sect. 13.2b . 
Each run starts from a snapshot of the detailed time-dependent 
(TD) simulation in the time window 600 s to 2350 s. The initial 
atmosphere of the individual runs contains the statistical equilib- 
rium (SE) solution for the thermodynamic state in that snapshot. 
In the next time step, the atmosphere is perturbed by increasing 
the gas temperature by 1 % and the time-evolution is followed 
for 50 s. In addition, we calculate the statistical equilibrium state 
of the perturbed atmosphere. For each snapshot in the dynamic 
simulation we thus have the following data: (i) the statistical 
equilibrium solution of that snapshot, (ii) the time-evolution af- 
ter a perturbation has been applied, and (iii) the final equilibrium 
state corresponding to the perturbed atmosphere. 

2.1. The calcium model atom 

The calcium model atom contains 5 bound states of singly 
ionised calcium (Call) and a continuum level (/ = 6), which 
represents the next ionisation stage (Calll). The lowest energy 
level with / = 1 is the ground state of Call (4s ^S). We also in- 
clude the two important energy level pairs of 4p^P (/ - 2,3) 
and 4d (/ = 4, 5). All five allowed radiative bound-bound (b- 
b) transitions are considered: the H and K resonance lines, and 
the infrared triplet (IRT). Also all five radiative bound-free (b- 
f) transitions with the corresponding photoionisation continua 
from the five lowest levels to level / = 6 are included. The model 
atom is illustrated in Fig.[T] 

3. Results 

3.1. Ionisation fraction 

In the solar chromosphere, calcium is mostly present in singly 
and doubly ionised form, i.e. as Call and Calll, respectively. 



2 



Sven Wedemeyer-Bohm and Mats Carlsson: Non-equilibrium calcium ionisation in the solar atmosphere 



-0.5 



-1.0 



-1.5 



-2.0 




O 



-2.5 



-4.0 



-5.1 (z„ = 1.4Mm) 



Ig m. 



-3.3 (z„ = 1.0Mm) 




: Ig m^= -1.3 (Zo = 0.5 Mm) 




1500 
time [s] 



Fig. 3. Logarithmic Call-III ionisation fraction as function of 
time at different Lagrangian locations. The labels refer to the 
corresponding geometric height in the initial atmosphere. The 
lines represent the results from the non-equilibrium simulation 
(solid) and for statistical equilibrium (dotted) calculated for the 
same values of the hydrodynamic variables. The horizontal lines 
are the corresponding time-averages. 



We therefore neglect neutral calcium and define the Call-Calll 
ionisation fraction as the ratio of the population density of dou- 
bly ionised calcium wcaiii and the population density of all in- 
cluded Ca ions (nca, total = "Caii + ncaiii)- For our six level atom 
it reduces to 



XC^L = 



ncaiii 



"6 



"Call + nCalll Xtl "i 



(1) 



where «, is the population density of level / and ncaiii - "e- 
In the following, we refer to xca as Ca II-III ionisation fraction. 
The variation of xca is shown in Fig. |2] as a function of col- 
umn mass density in the initial radiative equilibrium atmosphere 
(thick line). It has a minimum of 2 10 in the middle photo- 
sphere at Igmc = 0.0, which corresponds to zo ~ 180 km on the 
height scale of the initial model. Like the hydrogen ionisation 
fraction (dashed line) it rises with height but less steeply. At the 
base of the transition region, xc-d ~ 5 % is reached. The grey 
lines in Fig. |2] illustrate the temporal evolution of the Ca II-III 
ionisation fraction during the simulation. The thin dotted lines 
are the 5 % and 95 % percentiles, which enclose the typical data 
range of xca- While the temporal variation is small below the 
middle photosphere, the fluctuations cover usually two orders 
of magnitude in the layers above. This is caused by the upward 
propagating shock waves. 

In Fig. [3] the Ca II-III ionisation fraction xca is shown as 
function of time for three different Lagrangian locations, i.e. 



at fixed column mass densities. The corresponding geometri- 
cal heights z vary in time. For orientation, the labels in each 
panel specify the geometrical height zq in the initial atmosphere. 
Please note that the fluctuations of the ionisation degree appear 
much larger at a fixed geometrical height compared to a fixed 
column mass density because the atmospheric stratification and 
with it the height of the transition region is influcenced by the 
propagating shocks. 

The ionisation degree in the time-dependent non-equilibrium 
simulation (TD, solid line in Fig.O is compared to the statistical 
equilibrium values (SB, dotted line in Fig. O. Both are calcu- 
lated for the same atmospheric states. The differences between 
both cases are hardly discernible in the lower panel, which refers 
to the top of the photosphere. It shows that the ionisation frac- 
tion does not deviate considerably from the statistical equilib- 
rium values in the lower atmosphere. The deviation is somewhat 
larger but still moderate in the upper chromosphere at heights 
around z ~ 1400 km (upper panel in Fig.[3]l. At times of minimal 
ionisation, the time-dependent and statistical equilibrium results 
are very similar. At phases of higher ionisation in connection 
with passing shock fronts, the time-dependent simulation gives 
slightly smaller peak values, which lag behind the SE solution 
by 5 to 10 s. The time-average of the SE case is ~ 27 % larger 
than the TD result (see horizontal lines). 

In the middle chromosphere at heights around z ~ 1000 km 
(see middle panel of Fig. O, the SE ionisation fraction (dotted 
line) varies over two orders of magnitude on time scales of a few 
tens of seconds. In contrast, the detailed time-dependent sim- 
ulation produces a Ca II-III ionisation fraction that cannot fol- 
low the rapid changes (solid line). After a strong increase and 
a shock-induced peak, which is similar in both cases, the ion- 
isation fraction decays so slowly in the TD simulation that al- 
ready the next shock front has arrived before the equilibrium 
value can be attained. As we shall see in Sect 13.31 this be- 
haviour is caused by a relatively slow net recombination pro- 
cess. There is also a general time difference of the order of 5 to 
10 s between the occurrence of the ionisation peaks in the time- 
dependent simulation compared to the SE approach. At heights 
around z ~ 1000 km, the TD ionisation fraction nevertheless 
varies strongly between values of ~ 10"^ and ~ 10 Although 
the temporal evolution of xc-d is very diff'erent for the TD and 
SE case at those heights, the time-averages are relatively similar. 
The SE average is only ~ 14 % smaller than the TD result (see 
horizontal lines). 

3.2. Timescales 

In the following, we derive the relaxation timescale of the 
Ca II-III ionisation fraction (i) numerically from the tempo- 
ral evolution of an initially perturbed atmosphere (Sect. 13.2.1b 
and (ii) from the eigenvalue analysis of the rate matrices 
(Sect. 13^ . 

3.2. L Numerical relaxation timescale 

The simulation runs are used to determine a numerical timescale 
as function of time by analyzing each time step in the time 
range from 600 s to 2350 s individually (see Sect. |2]i. For each 
time step, the timescale is derived from the exponential relax- 
ation of the Ca II-III ionisation fraction towards the equilibrium 
state of the perturbed atmosphere. The time-averaged relaxation 
timescale is then calculated by averaging over the timescales of 
the individual time steps. This average timescale is shown in 
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Fig. 4. Relaxation timescale (solid lines) and gas tempera- 
ture (dot-dashed) as function of column mass and time in the 
dynamic simulation. The thick lines represent the relaxation 
timescale and the temperature averaged in time over all time 
steps between f - 600 s and f - 2350 s. Two individual time 
steps are drawn as thin lines: f = 1600 s (grey) and t = 1620 s 
(black). The dotted vertical lines mark the temperature peak and 
the minimum timescale for this time step. For reference, the 
height scale (zo) of the initial state is given on top. 



Fig. 5. Numerical relaxation timescale as function of column 
mass determined numerically (solid line) compared to the 
timescale determined from eigenvalues of the rate matrix in the 
initial statistical equilibrium states (dot-dashed). The thick lines 
represent the averages over all time steps between t = 600 s and 
t - 2350 s, whereas the thin lines are for a particular time step 
at f = 1620 s. Minimum and maximum values for both dynamic 
and eigenvalue timescales during the considered time window 
are shown as grey dotted lines. For reference, the height scale at 
t - 1620 s is given on top. 



Fig. |4] as function of column mass (thick solid line). It is very 
small in the photosphere (< Is) and increases strongly with 
height until a maximum with values of the order of 1 min is 
reached at heights around z = 900 km. The individual time steps 
show exactly the same behaviour in the photosphere, whereas the 
maximum in the chromosphere can reach values of up to 150 s 
in some cases. In the middle chromosphere above, the timescale 
is decreasing towards a minimum for all time steps. The average 
timescale drops to ~ 10 s at zo = 1200 km. The individual time 
steps also show a minimum but the exact position varies owing 
to the passage of shock waves. Values down to only a few sec- 
onds are found. The timescale starts to rise again at even larger 
heights but gets very small at the high-temperature transition re- 
gion. 

A number of successive individual time steps are represented 
by thin lines in Fig. H] They demonstrate the response of the re- 
laxation time scale (solid lines) to the passage of a shock wave 
(gas temperature as dot-dashed lines). At first glance, it seems as 
if the relaxation timescale and gas temperature are neatly anti- 
correlated, i.e., the shortest time scales are found at the posi- 
tions with the highest gas temperatures. This would be expected 
from high temperatures resulting in large ionisation rates and 
thus small timescales. 

A closer look, however, reveals that the minimum timescale 
is found at a height slightly below the temperature peak. This 
can be seen at f = 1620 s (thin black lines in Fig.lUi. The height 
difference is marked by dotted vertical lines. At the beginning of 
the shock wave passage in the photosphere, when it is still a dis- 
turbance with small temperature amplitude, the timescale at the 
position of the disturbance is mostly set by the thermodynamic 
state of the background atmosphere and thus the aftermath of 
previous shock waves. It is only at heights of z > 900km, where 



the shock has steepened enough to have a significant effect on 
the timescale. There, a local timescale minimum develops at the 
height of the temperature peak. While the shock front continues 
to propagate upwards in the atmosphere, the peak temperature 
grows, resulting in shorter and shorter timescales. The timescale 
minimum starts to appear lower than the peak temperature. This 
height difference grows during the passage through the chromo- 
sphere until values of the order of Az ~ 100km are reached. 
From the time at which the locations of the temperature peak 
and the timescale minimum coincide, the relaxation timescale 
drops from a few 10 s to just a few seconds. This behaviour is 
commonly found for all shock events. 

When looking at the data as function of time at a fixed ge- 
ometric height, the height offset translates into a small time dif- 
ference between the occurrence of the temperature peak and 
the minimum timescale. This time difference decreases in most 
cases with height from a few 10 s around z ~ 1000 km to only 
a few seconds around z ~ 1300km. There are examples where 
the time difference does not show such a height dependence but 
rather varies around values of just a few seconds. The time and 
height offsets are caused by the influence of the hydrogen ion- 
isation on the recombination rates of calcium via the electron 
density. See the analysis of the transition rates in Sect. 13.31 for 
more details. 

3.2.2. Timescales derived from eigenvalue analysis 

As in paper I, we use the eigenvalues of the transition rate ma- 
trix to determine the involved timescales and to indicate the pro- 
cesses governing the timescales. The rates in the dynamic sim- 
ulation are affected by the previous history of the atmospheric 
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state. The corresponding timescales are therefore not directly 
comparable to timescales derived from the numerical experi- 
ments in Sect. 13. 2. T] Therefore, we use the statistical equilibrium 
state corresponding to the hydrodynamic variables of a given 
snapshot (see Sect.|2]l. 

The Lagrangian time-derivative of the population densities 
is given by the rate equations: 



Dn 

— ^Pn 

Dt 



(2) 



with n being the vector of the level populations and P the rate 
matrix. The matrix element Pij of the rate matrix is the transition 
rate per atom from level / to level j with P„ = - 2; Pij- 

The matrix elements depend on the population densities in 
general through the non-linear radiation terms. If these non- 
linearities are not dominant we may write the solution to Eq. (|2]i 
in terms of eigenvec tors and eigenvalues of the rate matrix 
(Paper I; lJudgdl2005h . For the Ca model atom (Sect. |2TT]l with 
5 bound states and a continuum level (/ = 6, Ca III), the solution 
for the relaxation process can then be written as 



a Vi e 



Ail 



(3) 



Here, v,- is the ith eigenvector of the rate matrix P, Aj is the cor- 
responding eigenvalue, and c, is a coefficient, which depends 
on the initial conditions. The equilibrium state corresponds to 
the eigenvector with a zero eigenvalue. All other eigenvalues are 
negative and represent the evolution of the level populations to- 
ward their equilibrium values. The longest timescale, on which 
the system relaxes towards the equilibrium state, is then given 
by the inverse of the smallest (in absolute value) nonzero eigen- 
value. 

In Fig.|5] the timescales from the eigenvalue calculation are 
compared to the numerically determined relaxation timescales 
described in Sect. 13.2.11 The timescales averaged over all snap- 
shots between t - 600 s and t - 2350 s match closely 
throughout the whole atmosphere. The eigenvalue and numeri- 
cal timescales also match very well for the individual snapshots. 
This shows that in the case of calcium ionization, the timescales 
do not depend strongly on the non-linear parts of the rate ma- 
trix and we have a convenient way of calculating the ioniza- 
tion/recombination timescales from an eigenvalue analysis of the 
rate matrix. 



3.2.3. Timescales of individual processes 

We now investigate which processes dominate the relaxation 
timescale by modifying the entries of the rate matrix prior to 
determining the eigenvalue timescale. The resulting timescales 
have qualitatively the same run with height but differ in their ab- 
solute value. Leaving out the collisional rates does not change 
the overall timescale in the middle chromosphere because the 
collisional transitions occur on timescales that are up to 5 or- 
ders of magnitude larger compared to the overall timescale. Also 
setting the radiative b-b transition rates to zero has no notice- 
able effect on the timescale. Setting the rates Rk, and /?6i for 
the radiative b-f transition between the ground level (i - 1) and 
the continuum (/ = 6, Ca III) to zero changes the timescale only 
slightly. 

We now set all matrix elements to zero except for the en- 
tries of individual b-f transitions. A rate matrix with non-zero 
entries only for Ri g and R(, i results in timescales of up to a few 



thousand seconds in the middle chromosphere. It implies that 
this transition, which connects the most populated energy level 
(i - 1) with the ionisation continuum, does not govern the over- 
all timescale. It is the relaxation between level / = 3 and the con- 
tinuum, which produces the shortest timescales, followed by the 
transition between level / = 2 and the continuum. Both transi- 
tions have timescales of the order of a few 10 s. The ratio of the 
timescale with just a single transition and the timescale for the 
full rate matrix is about 2 for the transition 3-6 and ~ 3 for 2-6. 
The ratios for the transitions with lower level 4,5,1 are ~ 20, 10, 
and ~ 70, respectively. We conclude that the transitions between 
the metastable levels (2 and 3) and the continuum are the most 
important for setting the timescale for ionization/recombination. 
Including only these transitions result in a timescale that is only 
~ 18 % longer than the timescale derived from the full rate ma- 
trix. 



3.3. Transition rate analysis 

In contrast to the hydrogen case discussed in paper I, the shortest 
timescales do not coincide exactly with the temperature peaks in 
the shock fronts (see Sect. 13.2b . In this section, we show that 
this behaviour is caused by the dependence of the recombina- 
tion rates on the density of free electrons, which is influenced 
by the ionisation/recombination timescale of hydrogen. For the 
Ca model atom considered here, the recombination process in- 
volv es a collision b etween a free electron and a Ca III ion (see, 
e.g.. lMihalaslll978L p.l30f). Calcium is only a minor electron 
donor in the solar atmosphere, whereas already small changes in 
the hydrogen ionisation fraction can lead to significant fluctua- 
tions of the electron density. As shown in paper I, the hydrogen 
ionisation/recombination timescale is on average on the order of 
one to several hours in the chromosphere but can be strongly re- 
duced in hot shock fronts. The timescales can then get as short 
as 10 s to 20 s. The resulting delayed release of electrons has a 
direct effect on the recombination rates of Ca so that the related 
timescales reach their minimum only shortly after the occurrence 
of the temperature peak. 

In the following, the role of the individual ionisation and re- 
combination processes is illustrated in detail for a typical shock 
that propagates through the model chromosphere at f ~ 2200 s 
in the dynamic simulation. Most other shocks in the simulation 
show the same behaviour In Fig.|6l the gas temperature, Call- 
III ionisation fraction, and transition rates are shown as function 
of column mass density for the selected shock event. The four 
different time steps (rows from top to bottom) represent differ- 
ent stages of the upward propagating shock wave, which is most 
clearly seen as excess in gas temperature (dot-dashed line in the 
left column). The corresponding increase in ionisation fraction 
(solid line in the left column) is comparatively subtle with a peak 
(if discernible) slightly below the gas temperature peak. 

The normalised net rate P^jj for a radiative transition be- 
tween level / and level j is given by 



Pui = 



"Ca.total 



(4) 



with the level population densities «, and rtj and the radiative 
transition rates Rjj (upwards, / — > j) and Rjj (downwards, j — > /). 
The bound-free transition rates Pi-jd are shown in the right col- 
umn of Fig. |6] together with the gas temperature as reference. 
The net rates are thus given per Ca ion with positive val- 
ues indicating net photoionisation and negative values indicating 
net recombination. Collisional b-f transitions have much smaller 
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Fig. 6. Call-III ionisation fraction (left col- 
umn, solid line, left axis), gas temperature (left 
column, dot-dashed line, right axis), and nor- 
malised rates (right column) as a function of 
column mass for four time steps (rows from top 
to bottom) in the dynamic simulation. The rates 
for the radiative transitions from the bound lev- 
els 1-5 to the continuum (/ = 6, abbreviated 
"c") are divided by the total Ca number den- 
sities. Positive rates indicate ionisation, nega- 
tive rates recombination. The sum of all radia- 
tive b-f transitions is represented by the thick 
solid line. The gas temperature is repeated in 
the right column as reference (dot-dashed line). 
The height axis for the individual time step is 
shown at the top of each panel. Please note 
the different data ranges for the individual time 
steps. 



rates and can be neglected for this analysis. The sum of the ra- 
diative rates (thick solid) is close to zero in the lower parts of 
the model atmosphere (lg;7j(; > -1.0). There, the rates are in or 
close to equilibrium with the rates of the individual transitions 
all being small. The sum of the rates is also close to zero in the 
upper parts of the model atmosphere below the transition region 
(Ig/Wc ~ -5.6) that have not yet been disturbed by the shock 
wave (e.g., Igmc < -3.5 in the upper row). 

In the shock front, the sum of the rates is increased, lead- 
ing to net ionisation (positive values in the figure). The largest 
contribution is due to photoionisation from the meta-stable 
level i - 3, followed by photoionisation from level / = 2. As the 
lower levels of these transitions belong to the same term and are 
energetically close together (see Fig.lTJ, the transition rates show 
a very similar behaviour The photoionisation from the ground 
level a - 1) is seen at all heights below the transition region 
with a local maximum in the shock front, where it has a small 
but positive contribution to the net ionisation. The normalised 
rate reaches much larger relative values above the shock front 
close to the transition region, compensating the recombination in 
most of the other b-f transitions. These large values are a result 



of the previous shock wave and long relaxation timescales for 
this transition (see Sect. l3.2.3l l. The transitions from level / - 4 
and / - 5 have much smaller rates with recombination prevailing 
throughout the model atmosphere. They both behave very simi- 
larly because the lower levels are energetically close together. 

As the shock propagates through the chromosphere, the nor- 
malised photoionisation rates increase significantly. While the 
sum of the rates (thick solid line) is of the order of lO"'* s"^ at 
the shock front ?A z ~ 900km (top row of Fig. |6]l, values of 
~ 2.5 lO"-' s"' are found at z = 1 150km only 50 s later (bottom 
row). The second row from the bottom (f = 2250 s) shows the in- 
stant shortly before the shock wave reaches the transition region 
(seen as sharp jump in gas temperature around IgWc = -5.6). 
Once the shock reaches the transition region (20 s later, bot- 
tom row), the rates decrease rapidly. Except for a very small 
and smooth photoionisation from the ground level (/ =1), all 
other transitions and with it the sum of the rates show recombi- 
nation in almost the entire model atmosphere from just behind 
the shock front down to the low photosphere. It demonstrates 
that the timescales on which the ionisation equilibrium is re- 
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Fig. 7. Temporal evolution at a fixed geomet- 
ric height of z = 1200 km. Top: total Ca num- 
ber density nca.totai and relative level popula- 
tions Hi/ntotai- Middle: net radiative transition 
rates P^jj for the b-f transitions between the 
bound levels and the continuum together with 
the overall sum 2^=o ^n6- Bottom: the involved 
radiative rates Rjj and Rjj in comparison to 
the density of free electrons n^]. The Call-III 
ionisation fraction xca (grey dotted line), the 
gas temperature (grey dashed line), and the in- 
versed relaxation timescale (grey solid line) are 
shown in all panels as reference. The rates for 
the transitions from the bound levels 1 - 5 to 
the continuum (; = 6) are divided by the to- 
tal Ca number densities. Positive rates indicate 
ionisation, negative rates recombination. 



stored are relatively short - except for a comparatively sluggish 
photoionisation from the ground level. 

In Fig. 121 the time evolution of the shock event is analysed 
at a fixed geometric height of z = 1200km. The vertical dot- 
ted lines divide the time span into a pre-shock phase, a phase I 
(when the shock occurs) and a phase II (during the passage of the 
shock wake). As the processes are complex and the timescales 
are small, the Ca II-III ionisation fraction, the gas temperature, 
and the inverse relaxation timescale are shown in all panels for 
reference purposes (all scaled to their full data range). The gas 
temperature (grey dashed line) shows a clear shock signature 
with a rapid increase, a peak of Tgas - 8543 K at f = 2248 s 
and a slower decay in the shock wake afterwards. The Ca II- 
III ionisation fraction (grey dotted line) seems to react on the 
shock with some delay and also much smoother It peaks only at 
t - 2258 s with a value of 3.2%, i.e. 10 s after the occuiTence 
of the temperature peak. The relaxation timescale t reaches a 
minimum of Tmin = 7.6 s even later at f = 2262 s. The inverse 
timescale t"' is shown as solid grey line in the figure. While the 
change in phase I occurs rapidly (with a slope of t"' being be- 
tween those of the gas temperature and the ionisation fraction), 
the post-shock evolution in phase II is much slower. Then (e.g. 
att- 2300 s) the timescale r is still reduced while the gas tem- 
perature and the ionisation fraction are already relatively close 
to their pre-shock state again. 



The uppermost panel of Fig. |2]displays the total Ca number 
density nca.totai, i-C- the sum of levels 1 to 6 (thick solid line, 
right axis). As the element abundance of Ca is constant in the 
simulation, it behaves exactly like the gas density. It decreases 
before the aiTival of the shock front and rises rapidly again in the 
shock front and more gradually in the shock wake afterwards. 

The level population densities n,, i.e. the number densities of 
Ca atoms with their valence electron in a specific energy level /, 
are divided by the total Ca number density «ca,totai in the upper- 
most panel. It removes the fluctuations due to the change in gas 
density and reveals the relative distribution of the Ca atoms over 
the 6 included levels. At mid-chromospheric heights, usually 
95 % or more of all Ca II ions are in the ground level (i - 1) be- 
fore and after the shock passage (i.e., in the cool phases). During 
the passage of a shock wave, the ground level is depopulated but 
usually maintains at least 70 % of nca.totai and thus remains the 
most populated level at all times. The population densities n, of 
all other levels with / > 1 increase during the shock passage as 
Ca II atoms are excited from the ground level into a higher level. 
The increase of the population densities of levels 2 and 3 cor- 
relates with the change in gas temperature, while levels 4 and 
5 react with a delay of 5 s, shortly followed by the continuum 
level (/ = 6). The density of the latter, n^, increases at most to 
a few percent of nca.totai- Most Ca atoms therefore remain singly 
ionised, even during the shock passages. The Calll ions are ac- 
tually by a factor of more than 3 less abundant than the Ca II ions 
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Fig. 8. Emergent intensity in the Ca II K line: temporal evolution 
of the line core intensity (top row) and time-averaged spectrum 
(bottom row) at the solar limb (yu = 0.05, left column) and at 
disk-centre (// = 1.0, right column). Each panel contains the re- 
sults from the non-equilibrium simulation (/k,td, solid) and from 
the statistical equilibrium solution (/k,se, red thick dashed) to- 
gether with the corresponding difference (/k,td - ^k,se)/^k,se) 
between both (solid grey line at the top of each panel). The 
grey areas in the bottom panels represent the full data ranges 
of all intensity profiles during the simulation sequence between 
f = 600 s and f = 2350 s. 



with an electron in level 3 alone. Ionisation is thus only a "minor 
process" compared to the bound-bound transitions within the Ca 
atom. The level populations of levels 4 and 5 on the other hand 
are much smaller than the Ca III number density. Even if all ions 
in these levels would ionise, the resulting rise of the Ca III pop- 
ulation would be negligible. Levels 4 and 5 are nevertheless im- 
portant as upper levels of radiative b-b transitions. 

The net rates Pi ^ for the radiative bound-free transitions are 
shown as function of time in the middle panel of Fig. Q They 
are again given per Ca ion with positive values indicating net 
photoionisation and negative values indicating net recombina- 
tion (see Eq. Q). The sum of the net radiative rates 2;=i Pxjd 
(thick solid line) gives the change of the Ca III number density 
and thus the Ca II-III ionisation fraction. Contributions from col- 
lisional b-f transitions are again small here. The sum exhibits a 
net (photo-)ionisation in phase I, which peaks shortly (here 1 s) 
before the gas temperature, and rapidly turns into net recombi- 
nation for phase II. While the shock wake is passing through the 
analysed height point, the recombination becomes weaker and 
the rates approach zero again. The b-f transitions from levels 3 
and 2 contribute most (in that order) and show the same temporal 



behaviour as the sum. The net radiative bound-free rate for the 
ground level, f i ie, exhibits net photoionisation throughout the 
whole shock passage. The rate increases with time and reaches 
a first hardly discernible maximum together with or shortly be- 
fore the gas temperature. It roughly stays at the same level and 
smoothly declines in phase II. At that time, it is the only transi- 
tion with a positive net contribution. The bound-free transitions 
for levels 4 and 5 behave very much oppositely, producing a net 
recombination the whole time. The absolute values, however, 
stay small and even below the maximum recombination rates 
of Px,i(, and f r 36 that are reached during phase II. 

The peculiar behaviour of the net rates Px,ib and Pr,36 is 
caused by a temporal delay between the photoionisation {njRji) 
and the recombination part {iiiRjj, cf. Eq. dU). The photoionisa- 
tion terms increase first and produce the strong net photoionisa- 
tion in phase I, while the recombination part reaches a maximum 
only about 20 s later At that time both parts are equal, result- 
ing in a zero net rate. Afterwards the recombination part, which 
lags behind, stays larger than the ionisation part, causing the net 
recombination in phase II. The differences in the rates Rjj are 
balanced by the number densities of the lower levels so that the 
recombination part is only slightly smaller than the photoionisa- 
tion part. In contrast, the recombination into level / = 1 is much 
smaller than the corresponding photoionisation, which produces 
the smooth behaviour in time. 

Finally, the radiative bound-free transition rates /?, 6 for / = 
1,2,3 are shown in Fig. |7] The photoionisation rates /?26 and 
/?3 6, which are almost identical due to the similarity of their 
lower levels, exhibit no sharp response to the passage of the 
shock wave. In contrast, the recombination rates 7?6 3 and Rf,2 
increase strongly during the shock phase I. The rate /?6 3 is larger 
than R(, 2 but both only differ by a constant factor. The strong in- 
crease of these rates correlates very well with the inverse relax- 
ation time scale (grey solid line) and the electron density (solid 
black line). It demonstrates that the relaxation timescale is set 
by the recombination into the bound levels 3 and 2, which is di- 
rectly depending on the density of free electrons. With hydrogen 
being the main electron donor in the chromosphere, the tempo- 
ral behaviour of the recombination rates /?6 3 and R(, 2 is therefore 
similar to the H II number density. This leads to the conclusion 
that the Ca II-III ionisation in the quiet solar chromosphere is 
governed by the finite timescales on which hydrogen reacts on 
strong temperature changes. 

4. Discussion and conclusions 

The ionisation/recombination timescales, which are here derived 
from 1-D RADYN simulations, are too long for the assump- 
tion of an instantaneous ionisation equilibrium to be valid. On 
the other hand, the timescales are not long enough to warrant 
an assumption of a constant ionisation fraction. We find no- 
ticeable deviations from the ionisation equilibrium in the mid- 
dle model chromosphere but in general the Ca II-III ionisation 
fraction remains small. The error due to the often made simpli- 
fying assumption of statistical equilibrium is therefore negligi- 
ble for most applications. The effect is barely visible in the syn- 
thesized intensity for the diagnostically important spectral lines 
of Call, i.e., the H and K lines and the infrared triplet. This 
finding is illustrated for the CallK line in Fig. [8] The differ- 
ences in the emergent intensity between the statistical equilib- 
rium approach (/k,se) and the time-dependent non-equilibrium 
simulation (/k,td) are very small and hardly discernible in the 
temporal evolution of the intensity at a fixed wavelength posi- 
tion in the line core (top row) and even smaller (order of 0. 1 %) 
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in the average spectrum (bottom row). The relative difference 
(^K,TD - ^k,se)/^k,se reveals peaks with values of up to a few per- 
cent for wavelengths close to the line core. This is true both at 
disk-centre (right column) and close to the limb (left column). 
The differences are even smaller away from the line core and the 
emission peaks. Noticeable deviations occur only in connection 
with shock fronts and thus sharp intensity increases. There, the 
effect is caused by a small temporal offset in the evolution of 
the level populations in the TD case with respect to the SE case, 
like it is seen for the ionisation fraction in Fig. |3] It is safe to 
conclude that the time-dependent non-equilibrium treatment of 
the Ca II-III ionisation appears to be of minor importance for the 
lower atmosphere in quiet Sun regions. 

One restriction of the simulations is the use of only one spa- 
tial dimension, which could cause the shock wave profiles to be 
too extreme with potentially too high peak temperatures. In 3-D, 
the shock fronts are expected to be weaker owing to the larger 
number of degrees of freedom. The consequences of the devi- 
ations from the Ca II-III ionisation equilibrium would be even 
less important in that case. Furthermore, we neglect the effect 
of incident radiation from the corona on the photoionisation of 
Ca. Most important in that respect would be the Lyman alpha 
line at /I = 121.5 nm. However, the Lyman a photons have too 
little energy to ionise Call from the ground state (11.88eV or 
A - 104.4nm). The photoionisation edges for levels 2 and 3, 
on the other hand, both lie close to the Lyman alpha line. The 
incident radiation certainly matters most at large heights close 
to the transition region. There, however, photoionisation from 
the ground level dominates, which is obviously not influenced 
by Lyman alpha photons. We conclude that the RADYN simu- 
lations employed here are sufficiently realistic for an evaluation 
of the importance of non-equilibrium effects for the ionisation 
of calcium in the context of a strongly varying chromosphere in 
quiet Sun regions. 
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